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We study charge transfer between donor-acceptor molecules subject to a mirror symmetry con- 
straint in the presence of a dissipative environment. The symmetry requirement leads to the break- 
down of the standard single reaction coordinate description, and to a new charge transfer theory, 
in the limit of low temperature, based on two independent reaction coordinates of equal relevance. 
We discuss implications of these results to charge transfer in DNA. 
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Charge transfer between large organic molecules in 
aqueous solvent plays a key role in biochemical reactions, 
particularly those involving animal and plant metabolism 
jl], ^]. Compared to electronic transport in solid state 
materials, electron transfer in proteins and DNA is char- 
acterized by low values of the tunneling matrix elements 
and strong coupling between the electronic and nuclear 
degrees of freedom ||. Charge transfer mainly occurs 
when the nuclear coordinates happen to adopt a value 
for which the donor (D) and the acceptor (A) state en- 
ergies are nearly degenerate. 

In the conventional Born-Oppenheimer, or adiabatic, 
description, such a degeneracy region transforms to a 
saddle-point in the energy landscape separating the D 
and A states, and the transfer rate reduces to the classi- 
cal Kramers theory for activated chemical reactions in a 
dissipative medium ||, ||. However, because of the low 
value of the tunnel frequency Ao for biomolecular charge 
transfer, the adiabatic assumption is often not valid and, 
as a result, a quantum-mechanical description is required 
to determine the charge transfer rate. The interaction of 
the nuclear degrees of freedom with the solvent medium 
plays an important role in this respect because not only 
does coupling to a finite-temperature heat bath allow ac- 
cess to the degeneracy points by thermal activation, but 
friction between the nuclear degrees of freedom and the 
dissipative solvent - unavoidable due to the fluctuation- 
dissipation theorem - increases the time spent by the nu- 
clear degrees of freedom in the degeneracy region [Q, || . 

The analytical theory of DA charge transfer between 
biomolecules is based on the assumption that, among the 
many nuclear degrees of freedom involved, there is al- 
ways a single dominant one: the reaction coordinate Y 
that leads across the lowest saddle-point in the adiabatic 
energy landscape separating the D and A states. Charge 
transfer takes place when Y is within the Landau-Zcncr 
distance l LZ oc hA a /AF §, @ of the de generacy point, 
with AF the typical force level on the nuclear degrees 
of freedom at the degeneracy point pTTj ] . Cross-over to 
the adiabatic regime occurs when the time spent in the 
Landau Zener region tlz is large compared to the typical 
tunneling time Aq 1 . 

Single reaction DA theory is a basic tool not only for 



describing donor-acceptor charge transfer processes but 
also for describing chemical reactions involving macro- 
molecules in general. This single reaction coordinate the- 
ory is assumed to be at least qualitatively correct, though 
relaxation of the remaining nuclear degrees of freedom at 
the transition state may renormalize the parameters of 
the theory. It is the aim of the present letter to demon- 
strate that the single coordinate description fails in the 
presence of symmetry constraints. To be specific, con- 
sider charge transfer between two identical molecules. 
The electronic degree of freedom will be represented by 
Pauli spin matrices with o~ z = 1 denoting the electron in 
the D state and a z = —1 in the A state. The electron is 
coupled to the same two nuclear degrees of freedom of the 
D and A molecule Y\ , respectively, Yi . The Hamiltonian 
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Here, M is the effective mass and V the potential energy 
of the two nuclear degrees of freedom. Two collections 
{x ai } of harmonic oscillators represent the environmental 
degrees of freedom. They are coupled separately to the 
nuclear degrees of freedom, and generate a frictional drag 
on Y\ and Y2 with a friction constant fll2fl: 
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In the absence of any symmetry constraints, the lowest 
order coupling between the nuclear and electronic degrees 
of freedom is of the form (Yj.— Y2)a z . Treating (Y±— Y2) as 
the reaction coordinate, we can apply the standard single 
coordinate formalism. However, if we impose a mirror 
symmetry Y — > — Y, {x a — » — x a } then this term is 
forbidden. Expanding the potential energy to the lowest 
order in the lowest order in the nuclear coordinates under 
the symmetry constraint gives a Landau-Ginsburg type 
potential: 



2 




FIG. 1: Figure illustrating two different reaction paths in the 
Yi , Y2 plane of the nuclear coordinates. The locus of degener- 
acy points between the donor and the acceptor state energies 
are indicated by the bold lines Yi = iYi. The Landau-Zener 
region, within which charge-transfer reactions can take place, 
is shown as a dashed line. The 2D Landau-Zener length £ 
defines a region surrounding the origin where single-reaction 
coordinate theory fails. At the center of the 2D region one en- 
counters strong resonant tunneling within a distance a from 
the origin. The top trajectory, with a large impact param- 
eter q, crosses the two degeneracy lines separately. Charge- 
transfer at the crossing points allows transfer to the acceptor 
state \R I) (dashed arrows). The bottom trajectory, with 
a low impact parameter, enters both the 2D region and the 
coherent tunneling region. 
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Here, k is a spring constant, v describes the stabilizing 
effect of the lowest order anharmonic term, and A is the 
dimensionless coupling constant between electronic and 
nuclear degrees of freedom. 

We restrict ourselves to the case 1 < A < 2. In this 
regime, the donor state energy surface (f |V| f) nas two 
minima, at Y\ = ±y/k(X — l)/v, Y 2 = 0, correspond- 
ing to a charge deformed molecule in one of two alter- 
native mirror related structures. We denote the 'left' 
and 'right' donor structures respectively by | L f) and 
I R |). The acceptor state energy surface (J. \V\ j) has 
two corresponding minima but rotated over 7r/2 in the 
(Yi,Y"2) plane, as shown in Fig. 1. To locate the de- 
generacy points, we must solve for (| |V| f) = (I 1^1 I), 
which yields two degeneracy lines Y\ = ±Y 2 that cross at 
the origin of the (Yi, Y 2 ) plane. The lowest degeneracy 
point, the putative transition state, is the origin with an 
energy barrier AE = jk(X — l) 2 /v. 

To see why a single reaction coordinate description 



is not consistent in the presence of the symmetry con- 
straint, we show in Fig. 1 two possible trajectories in 
the (Yi, Y 2 ) plane for a thermally activated hop in the 
donor state from the L to the R configuration. The tra- 
jectories are classified by the impact parameter q, the 
distance of closest approach to the origin. A typical tra- 
jectory will cross each of the two degeneracy lines at least 
once. At each crossing point a charge transfer event can 
take place, with the system eventually ending in one of 
the two acceptor minimum states. Within a single reac- 
tion coordinate description for the charge transfer events, 
each degeneracy line is at the center of a Landau-Zener 
region of width Ilz(q) oc TiAa/kXq. For low impact pa- 
rameters, the two Landau-Zener regions of the degener- 
acy lines start to overlap when Ilz(q) drops below q, i.e. 
when q is less than £ = ^/hA Q /kX. As a result, there 
is a region of size £ surrounding the origin where charge 
transfer is inherently two dimensional. If the thermal en- 
ergy is low compared to the barrier height AE, then this 
2D region will dominate the charge transfer reaction and 
the single reaction coordinate assumption is not valid. 

In order to analyze a 2D charge transfer event wegen- 
eralize the Smoluchowski-Zusman (SZ) method |@, |l3) to 
two dimensions. In principle, this method is restricted 
to the case of strong damping, when the thermal energy 
fcgT is large compared to hu) c , with lo c the relaxation rate 
of the reaction coordinates. For the ID case though, the 
SZ method reproduces the weak coupling non-adiabatic 
result in the limit of low tunneling rates. 

In a 2D SZ description, the 2x2 density matrix rijj 
of the "spin" degrees of freedom (i,j now indicate spin 
up and down respectively) obeys the following transport 
equation: 
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with Fij — — V (i\V\j) the 2x2 force matrix derived 
from Eq.Q. The first term on the right hand side of 
Eq. (Q) is the Smoluchowski operator for the 2D diffu- 
sive motion of a classical particle in a force field. The 
second term describes the precession of a spin 1/2 de- 
gree of freedom with the spin Hamiltonian H a given by 
Eq. (||). Note that the energy difference between spin up 
and spin down has a saddle point at the origin. The 
characteristic length scale appearing in H a is the 2D 
Landau-Zener length £ encountered above. A compar- 
ison of H a with the Smoluchowski operator produces the 
second length scale a 2 n = {nksF '/rjk)i , which is the 
characteristic variation length of m,a- In the regime of 
strong damping, a 2 D is small compared to £. 

The equation for the off-diagonal part of the density 
matrix can be solved under the condition that n\ 2 varies 
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over length-scales that are short compared to those of the 
diagonal terms mi and n 2i2 : 

/mn li2 (?,t)~4^ 2 [n 1 ,i(0 1 t)-n 2a (0,t)}8(Y). (6) 

The 2D delta function in Eq. (|J) for the nuclear degrees 
of freedom is actually a gaussian with a width of order 
Inside this gaussian region, the off-diagonal part 
of the density matrix is large: the spin degree of free- 
dom is precessing coherently. Physically, this means that 
resonant tunneling is taking place between the D and 
A states within a distance ol 2 d of the origin of the (Y± , 
Y 2 ) plane. Substituting Eq.@ in Eq.(|), we find that 
the diagonal terms of the density matrix obey a classical 
Smoluchowski equation with a "sink" at the origin: 
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The decay rate of the donor state can be computed from 
Eq. (||) and Eq. (^|) using standard methods: 
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Here, w" 1 = rj/k(X — 1) is the classical life-time of the 
transition state, and To is the classical Kramers rate for 
an activated hop to the transition state. The constant C 
equals 4^2 In 2. 

The 2D decay rate of Eq. (^J) is in the same form of the 
ID description for donor-acceptor charge transfer pro- 
vided we interpret g 2B > — ftA^/ uj c k B T as the new adia- 
baticity parameter. For a ID charge transfer event, the 
adiabaticity parameter is of the form gn> cx fiA^ /lu c AEi, 
with AEb the characteristic energy scale of the nuclear 
degrees of freedom, such as the activation barrier. In 
the low temperature limit k B T « AE^, the effective 2D 
adiabaticity parameter diverges while gin remains finite. 
Physically, this means that in the low temperature limit, 
the adiabatic description is always valid in the presence 
of the symmetry constraint. 

We can justify treating g 2B) as an adiabaticity param- 
eter by estimating the time t 2B > spent in the 2D Landau- 
Zener region during a hop event. By Einstein's relation, 
the classical diffusion constant D of the nuclear degrees 
of freedom is k B T/rj, so that t 2B) ~ £ 2 /D is of order 
£ tj/^bT. Using £ = ^JhA^/k\ and identifying Aqt 2B > 



as the adiabaticity parameter we recover the above ex- 
pression for g 2 D- Interestingly, even though the quasi- 
classical method should become generally valid in the 
low temperature limit due to the symmetry constraint, 
the latter enhances the resonant tunneling regime as well. 
In the ID description, resonant tunneling takes place 
within a distance a\ B = {k B TTi/'qAF) 1 /' i of the degen- 
eracy point. Compared to a 2B , = (k B Th/ '^fc) 1 / 4 for the 
present case, we see that for T — > the regime of coherent 
tunneling in 2D is always larger than in ID. 

We conclude by noting that, a DA charge transfer sys- 
tem can be viewed as a two-state quantum system cou- 
pled to quasi classical macroscopic degrees of freedom 
(the nuclear coordinates) . Systems of this type form the 
basic elements for quantum computation, the "qubits" . 
For our result to be of any interest in this context, the 
adiabaticity parameter Aqt^z must be very large com- 
pared to one, of the order of 10 4 or larger [ fT5[ . Since the 
symmetry constraint leads to the divergence of the adia- 
baticity parameter in the low temperature limit, DA sys- 
tems obeying such a constraint would be better suitable 
for that purpose. A possible realization of the symmetry 
constraint is encountered for charge transfer between two 
DNA bases. Molecular dynamics simulations by Chen 
et al. [pr[ report that the in plane mirror symmetry of 
individual DNA bases is broken when a hole occupies 
the purine base. The mirror symmetry for stacked bases 
along a B DNA chain is not exact due to the propeller 
twist but the two alternative L and R structures of a 
charged base still may be sufficiently close for the present 
model to be relevant. It would be very interesting to test 
whether charge transfer events along DNA correlate with 
thermally activated hops between the L and R states. 
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